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Many models have recently been proposed in which dark matter (DM) couples to Standard Model 
fields via a GeV-scale dark sector. We consider scenarios of this type where the DM mass, at the 
electroweak/TeV scale, is generated by the VEV of a singlet which also couples to the Higgs. Such a 
setup results in a distinct recoil spectrum with both elastic and inelastic components. We construct 
an explicit NMSSM-like realization of this setup, discuss constraints coming from the relic density, 
and include benchmark points which are consistent with current limits, yet visible at upcoming 
direct detection experiments. 
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I. INTRODUCTION 



By now, there is overwhelming evidence supporting the 
existence of particle dark matter (DM) . Many models in- 
corporate DM into proposals for new TeV-scale physics, 
with recent efforts focusing on the possibility that DM 
may also interact with a dark sector, composed of new 
gauge groups and light (GeV-scale) degrees of freedom. 
These models attribute positron/electron cosmic ray ex- 
cesses [I] to the annihilation [2HS] or decay of [3 |H] of 
DM via the new GeV-scale states. In effect, the dark 
sector serves as a restrictive portal, allowing DM to de- 
cay/annihilate to light leptons, but not into (unobserved) 
protons and anti-protons. Dark sectors can also natu- 
rally induce an 0(100) keV mass splitting between the 
DM states [21 [3]. Consequently, if the fields mediating the 
scattering of DM with atomic nuclei couple off-diagonally 
to different DM states, the result is a novel inelastic recoil 
spectrum ^ which has been used to explain the DAMA 
results [To]. Indeed, regardless of any DAMA signal, in- 
elastic scattering is a generic consequence of dark sectors 
which can be probed by upcoming experiments. 

The phenomenology summarized above follows from 
the connection between GeV-scale fields and the Stan- 
dard Model (SM). But this is unlikely to be the whole 
story, because we would like the TeV-scale mass of DM to 
be related to the scale of electroweak symmetry breaking. 
The simplest way to generate the mass of DM is to cou- 
ple it to a singlet that also couples to the SM- Higgs and 
receives a TeV-scale VEV [2 [TT] . This naturally allows 
DM to scatter elastically off of atomic nuclei via the ex- 
change of a Higgs/singlet. We find that such dark sector 
models yield a distinct recoil spectrum, with both elastic 
and inelastic components, visible at the next generation 
of direct detection experiments such as XENONIOO. 

This paper is structured as follows. In Sec. [IT] we in- 
troduce this new recoil spectrum and survey its unique 
features. In Sec. |III| we construct an explicit NMSSM-like 
model realizing the scenario we propose, and discuss the 
constraints imposed on it from considerations of the relic 



density. In Sec. IV we discuss the masses and couplings 
of the model in a convenient limit, and present several 
benchmark points. Sec. |V] contains our conclusions. 
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FIG. 1: Sample semielastic recoil spectra with Xenon, corre- 
sponding to the two benchmark models of SectionllVj We have 
assumed a local DM density of p = 0.3 GeV/cm'^, and DM 
velocities are taken to follow a truncated Maxwell-Boltzmann 
distribution with Vq^^^ = 220 km/sec and a cutoff of Vesc = 600 
km/sec. We use the Helm nuclear form factor [12], with the 
parameterization of Ref. [T^ . The double hump structure on 
the left results from the inelastic component of the spectrum 
overlapping the zero of the form factor at ~ 100 keV. The 
summer (winter) spectra correspond to June 2nd (December 
2nd). Both spectra are consistent with current limits jl4[ I15j 
and predict ~ 20 events in the first year of XENONIOO data 
(within 8.7 o 40 keV). 



II. SEMI-ELASTIC SCATTERING 

The elastic scattering of DM with atomic nuclei can be 
a natural consequence of the mechanism setting its TeV- 
scale mass. We will illustrate this in an especially simple 
setup by coupling DM and the Higgs to the same singlet 
field. While we will work in a supersymmetric framework, 
our results are easily generalized to other scenarios. 

Consider the superpotential 



W = \SHi ■ Hu + vSxX, 



(1) 



where Hd and iJ„ are the two Higgs doublet fields, S is 
the NMSSM singlet [16 , and x/x, which are oppositely 
charged under a dark sector gauge symmetry, will com- 
pose our DM candidate. After electroweak symmetry 
breaking, S receives a VEV, which generates a super- 
symmetric mass for the DM fields, — r]{S). Taking 
DM to be a scalar component of x/x will justify this 
assumption in the next section), we see that the F-term 
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potential includes a direct coupling between DM and the 
Higgs, D xxHuHl + h.c. This coupling allows the 

Higgs to mediate the elastic scattering of DM against 
nuclei. A second contribution to elastic scattering is me- 
diated by the singlet S*, which mixes with the Higgs after 
electroweak symmetry breaking. 

At the same time, other dark sector interactions can 
lead to an inelastic component of scattering. We con- 
sider the class of models where DM is charged under a 
GeV scale dark sector, with a U{\)d gauge factor that 
kinetically mixes with hypercharge, 

^ ^ ^B^'^b^. + 9d b^'ixod.Xi - Xid^Xo)- (2) 

Here xo,i are the two real scalar components of DM sep- 
arated by a small mass splitting of order S ^ 100 keV, i?^ 
and 6^ are the hypercharge and U{l)d gauge fields, and 
e parameterizes the size of the kinetic mixing. The mass 
splitting can be generated by a higher dimension opera- 
tor [5], or radiatively by the breaking of a non-Abelian 
dark sector gauge symmetry [S]. Through kinetic mix- 
ing, the dark sector photon, b^, acquires e-suppressed 
couplings to quarks and thereby mediates inelastic scat- 
tering between DM and nuclei. This scattering will take 
place along with the elastic scattering described before, 
realizing a scenario we term semielastic scattering |27j . 

Phenomenologically, semielastic DM can be parame- 
terized by the 4 parameters, m^^, S, <te, and cr/. The 
elastic scattering [ue) dominates at low nuclear recoil 
energy while the inelastic scattering (cr/) dominates at 
higher recoil energy. Examples of such a spectrum are 
shown in Fig. [l] The unique spectral shape changes the 
constraints and reach, in the (cr£;,m;^o) plane, as shown 
in Fig. [2| We do not attempt to fit the possible DAMA 
signal [TU] . 



III. MODEL 

We now provide an explicit realization of the scenario 
described above. We take as our starting point the 
NMSSM, where a singlet superfield S couples to the two 
Higgs multiplets of the MSSM. To this, we add an addi- 
tional coupling of the singlet to DM, as in Eq. [T] and a 
U(l)rf dark sector along the lines of Refs. [51 ITS]. 

In detail, we will be concerned with the following terms 
in the superpotential, 

W D \SHd- Hu + 1lSxx+\^^S^ + pNRR+\x^R' (3) 

where Hd and Hu are the two Higgs doublet fields, S is 
the NMSSM singlet, x and x '^iU compose our DM can- 
didate, R and R are GeV-scale dark sector Higgs fields, 
and iV is a GeV-scale singlet whose presence insures that 
all dark sector fields receive a tree- level mass |6] . We will 
see below that the higher-dimension operator, suppressed 
by A > 10 TeV, will generate a small DM mass splitting. 



Semielastic DM Inelastic DM 




DM Mass [Gt:Vl DM Mass [GcV] 

FIG. 2: The left panel shows the current 90% limits, from 
CDMS [H] and XENONIOO [15], on the semielastic scenario 
in the DM mass - elastic cross-section plane. We have fixed 
the DM splitting, 5 = 140 keV, and inelastic cross-section, 
aj = 1.7 X 10"^" cm'^ per proton, to match the second bench- 
mark of Section [TVl We have also included the projected limit 
from XENONIOO after one year of data, assuming zero back- 
ground, a raw exposure of 6000 kgxdays, an efficiency of 38%, 
and a nuclear recoil energy range of 8.7 to 40 keV. The promi- 
nent dip in cross-section indicates the range of masses where 
the model will be visible from inelastic scattering alone (the 
dashed orange curve shows the limit without inelastic scat- 
tering). This can be seen on the right panel which shows the 
current and projected limits on inelastic scattering only, with 
5 — 140 keV. The black (dashed) horizontal line on the right 
corresponds to the inelastic cross-section assumed on the left. 

We assign x and R (x and R ) charge 1 (-1) under a 
dark V{l)d gauge group. Furthermore, we assume the 
presence of supersymmetric kinetic mixing, 

cd-^J d^ewyWb, (4) 

for Wy and Wb, the hypercharge and dark supersym- 
metric field strengths, where e ~ 10^'' ^ 10^^ is natu- 
rally generated, at one loop, by integrating out physics 
at higher energy scales. Expanding in components, the 
kinetic mixing includes D-term mixing, which generates 
an effective Fayet-Iliopoulos D-term in the hidden sector 
at the GeV scale [Hl^. 

Upon minimizing the dark sector potential one finds 
that Rc develops a VEV (Re) = Vr ^ GeV which Higgses 
the dark photon, giving it a mass m^^ ^ dd^r- The other 
light dark sector states also live at the GeV-scale. Finally, 
we note that there can be 0(1) corrections to their masses 
coming from SM SUSY breaking, which is communicated 
to the dark sector through gauge interactions with x act- 
ing as a messenger. These corrections have been included 
in our benchmark spectra of section |IV[ although they 
have no qualitative effect on the phenomenology. 

We now consider the spectrum of the x rnultiplet, 
which will contain DM. After S gets a VEV, there 
are two nearly degenerate fermionic states with masses 
^ rjVs/V^ (these states are split a small amount by the 
higher-dimension operator). Meanwhile, the scalar com- 
ponents are split from their supersymmetric masses by 
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(Fs) , and under the assumption that x's dominant source 
of SUSY breaking is communicated by S", we can neglect 
additional soft terms. The four scalar degrees of freedom 
divide into two pairs with masses above and below the 
fermions, separated by a large weak-scale splitting, 



m = T] 



(77±k)=f 7w|w sin(2/3) 



(5) 



where Vs is the singlet VEV. 
there is a smaller splitting 



Within each scalar pair 



VsVy 

A 



(6) 



where Vr is the VEV of the dark Higgs. In what follows, 
we will label the scalar mass eigenstates Xi for i : — >■ 3 in 
order of ascending mass. The lightest state, Xo; will serve 
as our DM candidate, and 8m — m^^ — rti^^ ^ 100 keV 
is a consequence of Eq. [6j 

Now, demanding that this model reproduce the ob- 
served relic abundance of dark matter places strong con- 
straints on the different couplings and VEVs. DM can 
annihilate via three competitive channels: (1) to Higgses 
and singlets, (2) via dark sector gauge interactions, and 
(3) to dark sector Higgses through the higher-dimension 
contact interaction of Eq. |3] DM has the observed relic 
density if these channels sum to have the correct annihi- 
lation rate for a thermal relic, 



2.5 X 10--* GeV" 



(7) 



We derive constraints on the model by demanding that 
no individual channel exceed this rate. To derive these 
limits, it is sufficient to consider self-annihilations of Xo- 
In general, there is also co-annihilation of xo with xi 
(and the x fermions for small enough (Fs)), but the fol- 
lowing estimates still apply up to 0(1) corrections. We 
do not consider constraints on Sommerfeld enhanced an- 
nihilations from cosmology [19] . Our limits are therefore 
conservative and also applicable to decaying DM mod- 
els 0. 

First we consider DM annihilations into the SM Hig- 
gses. Assuming these are light compared to m^, one finds 
their contribution to the annihilation rate to be 



> 1 f 5"^xo 



8tt \ V 



A4 
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which implies 



w > 4.5 TeV 



VI TeV/ 



1/2 



(8) 



(9) 



While there are also constraints on A, they are far less 
severe: 



A < 0.4 



VI TeV/ 



1/2 



(10) 



DM can annihilate through dark sector gauge interac- 
tions with rate. 



9i 



Xo 



which constrains the dark gauge coupling: 



' 9d 
.0.5 



2 / 1 TeV 



< 1. 



(11) 



(12) 



In addition, the non-renormalizable operator that gen- 
erates the small DM splitting allows DM to annihilate 
into pairs of dark-Higgses, 



(av) 



(13) 



Therefore, we find a non-trivial constraint relating the 
mass splitting between the dark matter states and the 
dark sector breaking scale [8|, 



100 keV 



1 Gcv y 

Vr J 



< 1 



(14) 



We conclude this section with a brief discussion of 
other important constraints on this model. We note 
that the excited state xi is long-lived and has a relic 
density that is constrained by inelastic down-scattering, 
but its density can be depleted in several ways as dis- 
cussed by Ref. [5D]. There are a number of additional 
constraints if one attempts to explain the cosmic ray 
anomalies through Sommerfeld enhanced annihilations as 
in Refs. [21 Sj. Dark matter can annihilate into hidden 
sector gauginos which decay to a dangerous amount of 
SM photons as discussed by fS^. This constraint is alle- 
viated if mgravitino ^ 1 GcV, Or by considering a more 
elaborate dark sector. There is tension from astrophysi- 
cal limits on neutrinos and photons from final state radia- 
tion, coming from the galactic center JJ_J. Alternatively, 
these astrophysical tensions are alleviated if the cosmic 
rays are produced by DM decays into the dark sector [S] . 



IV. BENCHMARKS 

It is convenient to consider this model in the various 
analytically tractable limits of the NMSSM. One finds, 
however, that whether one starts with a small k (the 
PQ-symmetric limit [55]) or with small A-terms (the R- 
symmetric limit [23] ) . the requirement of a sizable Vs, sl 
stable EWSB minima, and an elastic recoil spectra visible 
at current direct detection experiments necessitates small 
K and A. We therefore consider the limit k, A — > 0. 

DM scatters elastically by exchanging the three CP- 
even Higgses, s, h, and H. With DM at the TeV scale 
and the three scalar Higgs masses above 100 GeV, we 
find that it will be difficult to see the elastic scattering 
at current direct detection experiments. Things become 
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TABLE I: Benchmark points for semielastic scattering. Here nis is the mass of the hghtest CP-even scalar, which tends to 
dominate the elastic scattering rate. The first benchmark has a DM mass that can explain the cosmic ray anomalies with 
annihilations (m^^ ~ 1 TeV), while the second benchmark has a mass appropriate for decaying DM (m^^p ~ 2 TeV) [21j . 
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more interesting if one of these states is light, enhanc- 
ing the elastic cross-section. The mostly-singlet scalar, 
s, enjoys suppressed couplings to the electroweak gauge 
bosons and can be very light without conflicting with 
existing LEP limits [23]. We work in the limit where 
rus ^ 50 GeV. One finds 

aei - 1.7-10-^° cni^ [ '"""^,''''"'' )' (15) 

where f;^ ~ 1 and qh ^ ^(tan/3 ~ cot/3) parameterize 
the couplings of h and H to nucleons [25, (we use the 
nuclear matrix elements of Ref. |26J), and ah/cnH denote 
the singlet-Higgs mixing angles. Note that while it seems 
one can arbitrarily increase aE via gn by choosing a large 
tan /3, the singlet proportionally decouples from i?, so no 
such tuning is possible. 

We present 2 benchmark points in Table |l] Both yield 
the correct relic abundance and a recoil spectrum visible 
in one year of XENONIOO data. 

V. DISCUSSION 

Here we have studied a mechanism that naturally re- 
lates the mass of DM to the scale of electroweak symme- 
try breaking in models with a GeV-scale dark sector. DM 
scatters against nuclei elastically via a Higgs/singlet, and 
inelastically through dark gauge boson exchange. Com- 



bined, the spectrum has a unique semielastic shape which 
can be discovered in upcoming direct detection experi- 
ments such as XENONIOO. 

While our primary concern has been the unique recoil 
spectrum particular to this class of models, we found in 
Sec. |III| that the parameters are constrained, nontrivially, 
by the requirement of getting the right relic density. It 
would be interesting to further investigate the interplay 
between these constraints and the modified Higgs phe- 
nomenology of the NMSSM. 
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